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Abstract. The entanglement content of high-dimensional random pure states is 
almost maximal, nevertheless, we show that, due to the complexity of such states, 
the detection of their entanglement using witness operators is rather difficult. We 
discuss the case of unknown random states, and the case of known random states 
for which we can optimize the entanglement witness. Moreover, we show that coarse 
graining, modeled by considering mixtures of m random states instead of pure ones, 
leads to a decay in the entanglement detection probability exponential with m. Our 
results also allow to explain the emergence of classicality in coarse grained quantum 
chaotic dynamics. 
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1. Introduction 

Random pure quantum states, i. e., states whose distribution is invariant under an 
arbitrary unitary transformation, are almost maximally entangled. Let us consider, for 
instance, a random state from N x N dimensional Hilbert space n = . We 
consider a bipartition to subsystems A and B, namely as an element of a tensor 
product Hilbert space Ha ® Hb where Ha = C^, Hb = are Hilbert spaces of 
subsystems A and B, respectively. The entanglement E of the state is then measured 
by the von Neumann entropy S{p) = — tr(plog2p) of the reduced density matrices 
Pa = trB|V')(V'| or pb = ti A\i'){'ip\- It turns out [HEIE] that 
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which is close to the maximal value log2 attained for maximally entangled state. 
Since random states carry a lot of entanglement and entanglement [H |5] has no analog 
in classical mechanics, one can conclude that random states are highly non-classical. 

On the other hand, pseudo-random states with properties close to those of true 
random states can be efficiently generated by dynamical systems (maps) in the regime 
of quantum chaos O [Tj El 13 [lOl [H] ■ In such chaotic maps the classical limit is recovered 
when oo. Therefore, one can argue that for high- dimensional random states, i.e., 
in the limit — oo, the quantum expectation value of an operator with a well defined 
classical limit will be close to its classical microcanonical average. According to this 
picture random states in a way "mimic" classical microcanonical density. Expectation 
values are therefore close to the classical ones. 

At first sight this is in striking contrast with the almost maximal entanglement of 
such states. How can we reconcile this apparent contradiction? In the present paper 
we are going to tackle this question by considering how can we detect entanglement 
of random states. By studying the detection of entanglement with decomposable 
entanglement witnesses we are going to argue that in the limit of large systems the 
detection of entanglement in a random state becomes increasingly difficult as it would 
demand the control of very finely interwoven degrees of freedom and a measurement 
resolution inversely proportional to the size of the Hilbert space, which seems hardly 
feasible experimentally. Therefore, as far as the detection of entanglement is concerned, 
high dimensional random states are effectively classical. 

Moreover, coarse graining naturally appears. For instance, one could repeat several 
times the measurement of the entanglement witness for a random state and the prepared 
random state would be different from time to time due to unavoidable experimental 
imperfections. We model this problem by considering mixtures of m pure random states, 
namely 



where the are mutually independent random pure states, but in general they are 
not orthogonal. We are going to show that the detection of entanglement is even 
more difficult for these mixed states, as it requires a number of measurements growing 
exponentially with m. 

There are other physical contexts in which formally the same kind of coarse graining 
(mixing of the state) naturally appears: (i) Time averaging. For example, if a state {ip) 
undergoes a time evolution \ip{t)) = U{t)\ilj) given in terms of some unitary dynamics 
U (t) , then the time average of a physical observable A over an interval T is essentially 
determined by expectation value tr (Ap) in the mixed state 



which has an effective rank m ^ T/tcorr where tcorr is a dynamical correlation time of 
the dynamics U(t). For a quantum chaotic evolution U(t), the state \ip{t)) can be, after 
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some time, arguably well described by a random state and the correlation time tcorr is 
expected to be short, so p ([3]) may be considered as a mixture ([2]) of m uncorrelated 
random states, (ii) Phase space averaging. Sometimes it is useful to represent quantum 
states in terms of distribution functions in the classical phase space, like the Husimi 
function (see e.g. [7]), which can be understood as a convolution of the Wigner function 
or its coarse graining over a phase space volume {27rhY in d degreees of freedom. In 
fact, the Husimi function of a pure state can be understood as a Wigner function of the 
following mixed state 

p = {271^,)-' j dgdpexp (-^K' + «~y)) T{q,p)\ij){ij\T^{q,p) (4) 

where T{q,p) are unitary phase space translation operators, and a is an arbitrary 
squeezing parameter. A random pure state has a Wigner function with random sub- 
Planck structures with phase space correlaton length /corr ~ ^ which is semi-classically 
smaller than the coarse-graining width ~ h}^'^, so p can be again considered as a 
mixture ([2]) of m random pure states with m ~ fr^/'^. 

The paper is organized as follows. In Sec. [H we review known results about 
entanglement witnesses and pure states entanglement. Depending on our a priori 
knowledge of the random state, two relevant cases can be distinguished: (i) We do 
not know in advance of which random state we are going to detect the entanglement, 
i.e., the random state is unknown. In such a case the best one can do is to choose 
some generic entanglement witness in advance, independently of the state we measure. 
Such situation can also be thought to arise in the case when we are not able to prepare 
an arbitrary entanglement witness but just some subset of witnesses, (ii) We know 
the random state in advance and we are able to prepare an arbitrary entanglement 
witness. In such case we can use the optimal entanglement witness for each random 
state separately. These two cases are discusses in Sees. |3] and HI respectively. Finally, 
in Sec. |5l we provide a brief discussion of our results. 



2. Entanglement witnesses 

First of all, let us introduce some known facts about entanglement witnesses and 
entanglement of pure states that we will need for the derivation of our results in 
subsequent sections. Given a pure state \ip), its bipartite entanglement content is 
completely specified by the Schmidt decomposition: 

= ^Pi\ai)(^\hi), (5) 

i 

where the Schmidt coefficients pi are positive real numbers satisfying Ylii /^f = 1 ^i^d 
\ai) and \bi) are orthonormal states on subspaces T^a and Tie- The squares of the 
Schmidt coefficients are equal to the eigenvalues of the reduced density matrices pa 
and pb. For such bipartition there are at most N nonzero Schmidt coefficients pi. 
For a state having r roughly equal nonzero Schmidt coefficients {pi ~ l/-\/r), the 
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reduced von Neumann entropy is S{pa) ~ log2 r. For instance, the simple (GHZ) state 
\tp) = -^lO . . . 0) + -^|1 . . . 1) leads to two nonzero eigenvalues of pa {pi = P2 = Vv^)' 
and simple Schmidt vectors |ai) ® \bi) = |0 . . .0) and \a2) ® I&2) = |1 • • • 1), resulting 
in S* = 1. On the other hand, a random pure state has all N eigenvalues nonzero, 
very complicated (random) Schmidt vectors |aj) and \bi), and almost maximal entropy 
S ~ log2 N — 1/ logg 4 in [2], [3]. Note, however, that all the eigenvalues of the reduced 
density matrix for a bipartite random state are very small. In fact, they decrease with 
A^, the largest being, on average, A/N while the smallest is [12j. The guiding idea 

of this paper is that this smallness of eigenvalues and the complexity of eigenvectors 
should be somehow reflected in the difficulty of detecting entanglement in such a state, 
despite the fact that the entanglement content of a random state is large. 

Besides pure states, we will also be interested in the entanglement of a mixture of 
pure states ([2]). A quantum state p of a bipartite system is called separable if it can be 
written as 

P =^PkPAk® PBk, Pk>0, ^Pfc = l, (6) 
k k 

where pAk and pBk are density matrices for the two subsystems. A state is entangled if it 
is not separable. To decide whether a given mixed state is entangled or not is a difficult 
problem [H [5] . Fortunately, there is an operational criterion which is able to detect the 
most useful entangled states. This is the famous positive partial transposition (PPT) 
criterion [13]: since separable states have positive semidefinite partial transpos^, all 
non-PPT states are entangled. Note, however, that for dimensions higher than 2x2 
and 2x3 there exist PPT-entangled states, known as bound entangled states [14] . 

A convenient way to detect entanglement is to use the so-called entanglement 
witnesses [H| [T5] . By definition, an entanglement witness is a Hermitian operator 
W such that tr (H^Psep) > for all separable states psep while there exists at least one 
state Pent such that tr (VTpent) < 0. Therefore, the negative expectation value of W 
is a signature of entanglement and the state pent is said to be detected as entangled 
by the witness W. The concept of entanglement witness is close to experimental 
implementations and detection of entanglement by means of entanglement witnesses has 
been realized in several experiments [El [171 [18] . Because it is easier to measure a larger 
negative expectation value of an entanglement witness, one can argue that the detection 
of entanglement is easier the larger this negative value is. The expectation value of 
W also provides lower bounds to various entanglement measures [191 120]. Estimation 
of entanglement entropy by measurement of observables has been considered in |21j . 

I Introducing and orthonormal basis {\ia) = \i)A <8) |q:)b} in the Hilbert space associated with the 
bipartite system, the density matrix p has matrix elements pia.jp = {ia\p\j(3). The partial transpose 
is constructed by taking the transpose of only Latin or only Greek indexes (here Latin indexes refer to 
subsystem A and Greek indexes to subsystem B). For instance, the partial transpose with respect to 
subsystem B is given by 

Pia,j0 ^ PiP^ja- (7) 
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Finally, it is interesting to note that violation of Bell inequalities can be rewritten in 
terms of nonoptimal entanglement witnesses [22] . 

In general, classification of entanglement witnesses is a hard problem. However, 
much simpler is the issue with the so-called decomposable entanglement witnesses (D- 
EW) [23J. D-EW is a witness which can be decomposed as 

W = P + Q^^, P,g>0, (8) 

that is, with positive semidefinite operators P, Q. D-EW can detect only non-PPT 
entangled states, i.e., those with negative eigenvalues of p^^. They are therefore 
equivalent to the PPT criterion but closer to experimental implementations, as full 
tomographic knowledge about the state is not needed. In the present paper we are 
going to limit ourselves only to D-EW. General non-decomposable entanglement witness 
(ND-EW) can be written in a canonical form as W = P + Q^^ — el [23l [2l] and can 
also detect entangled states with PPT. Finding an optimal ND-EW, for which violation 
of positivity is maximal, is in general hard [23] . 

3. Unknown random state 

In this section we assume that the random state {ip) whose entanglement we would like 
to detect is unknown so that we are not able to use an optimal W for a particular 
The best one can do is to choose some fixed entanglement witness W in advance, 
independently of the state. Since we will be interested in the average behaviour over 
unitary invariant ensemble of pure random states, W can be chosen to be random as 
well. That is, in the present section we are going to study detection of entanglement 
with a random entanglement witness, whose precise definition will be given later. What 
we want to calculate is the distribution of the expectation values for a fixed 

W and an ensemble of random pure states {ip). Averaging over random states we 
see that the average expectation value (t/'I is 

{^\W\^) = j dV{ip\W\^) = tiW/N^, (9) 

where • = / dV» denotes an integration over a unique ?7(A^^)-invariant distribution of 
pure states lip), and we used the fact that for a random state \ip) = Xli'^iK) have 
CjC* = Sij/N"^. We fix normalization of the entanglement witness W such that tiW = 1. 
Therefore, the average expectation value (t/'I VTl?/') scales oc l/N"^. We therefore define 
the rescaled quantity w = N'^{4'\W\ip), such that w = 1, independently of the dimension 
A^. From now on we will focus on the random variable w, and its distribution with 
density p{w) = dV/dw. 

Because operator P in D-EW (IHl) just shifts the expectation value towards positive 
values, we limit ourselves to D-EW of the form W = Q"^^. Any positive semidefinite 
operator Q can be written in its eigenbasis in terms of positive eigenvalues di, satisfying 
— 1' eigenvectors hence W = ■ ^^^^ ^^^t study 

the case when Q is a simple rank one projector (subsection 13. ip . that is W is given by 
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W = There are two reasons why this is the most important case. First, 

as we will see below, optimal D-EW is always of such a "projector" form. Second, the 
expectation value for a general Q can be written (subsection [32]) as a sum of expectation 
values for individual eigenvectors of Q and therefore the probability distribution of w 
will be a simple convolution of distributions for the case of Q being rank one. 

Most of our theoretical results are derived for one mixing component only, m = 1 
(i.e. for a pure random state), since the general result for arbitrary number m of 
independent mixing components is obtained by simple convolutions as discussed in 
subsection 13.31 

3.1. Q is a projector 

The entanglement witness is of the form W = (|0)(0|)'^^, with a fixed and we would 
like to calculate the density of probability distribution p{w) of its normalized expectation 
values w for random pure states. This distribution can depend on the chosen First 
few moments of the density p{w) can be calculated explicitly. We have already seen that 
W = 1. For higher moments we can in the leading order (in Hilbert space dimension 
l/N"^) use Gaussian averages and Wick contractions in order to approximate integrals 
over the unitary group [25]. Using 

= {5,,5ui + 5,i5,k)/N^ + 0{N~^), (10) 

and similarly for higher order products, we arrive, up to corrections (9(A^^^), at 

K2 = {w - wy = trW^ = 1, 

= (w-wf = 2tr^y^ (11) 

= {w- wy - 3kI = 6tTW\ 

Here we denote by /t„ the n-th cumulant. We see that the average value of w as well as 
the width of the distribution is 1, independently of the state |0) we use for W. While 
tr FT = 1 is a simple normalization, second moment is 1 due to Q being of rank 1. If only 
the first two cumulants were nonzero our probability density p{w) would be a simple 
Gaussian. 

Let us first see what happens if we choose for |0) a state with large Schmidt rank 
r (number of nonzero eigenvalues of the reduced density matrix a a = trB|0)(0|) of 
order r ~ A^. Using Schmidt decomposition ([5]) for the state |0) = ® 
we can immediately write eigenvalues and eigenvectors of the Hermitian operator 
W = (10) (01)'^''. There are A^(A^ — 1) eigenvalues ±fiifij, i < j, with the corresponding 
eigenvectors (|aj6*) ± \ajb*))/\^ and eigenvalues fij with the corresponding 
eigenvectors \aib*). In our notation, \aib*) = \ai) ® \b*), where \b*) = ^Q&*Q,|a) if 
\bi) = ^Q,&ia|a). Using these eigenvalues one can see that the traces of powers of W 
are 

trW''^ = {tiaXf, trW^^+^ = tiaf+\ (12) 

Because eigenvalues of a a are of order ~ 1/r, we have trPi^^ ~ 1/r^ ~ (9(1 /A^^), 
triy^ ~ 1/r^ ~ (9(1/A^^) and so on. Note that formulas (ITTl) are exact up to order 
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and therefore can not be used for k„>3 and states |0) with the full Schmidt rank. 
Nevertheless, we can conclude that the higher cumulants are at most ~ l/N"^ and 
therefore vanish in the limit N ^ oo. Therefore, for W and |0) with an increasing rank, 
the probability density p{w) converges to a Gaussian in the limit — > oo, 

1 . . 



p{w) 



exp 



[w 



1)72). 



(13) 



Theoretical prediction f|T3|) is compared with the results of numerical simulation in 



figure [H Probabihty to measure negative w, 
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Figure 1. Density of probability distribution of w — N'^ {ip\W\ip) for random vectors 
and a single W = (|0) (</'|)'^^ with a random Total Hilbert space dimension 
— 2^". Full curve is the theoretical Gaussian prediction pS]). We use an ensemble 

of 5 X 10^ random states. 

J^^p{w)dw and therefore 

V(w < 0) = (1 - erf(l/V2))/2 ^ 0.159. (14) 

Note that this entanglement detection probability is independent of the details of 
provided that its Schmidt rank r is large, more precisely r oc N. 

The above Gaussian form of p{w) can be understood also as a consequence of the 
central limit theorem. Indeed, let us write W in its eigenbasis, then we have, writing 

■^i = /^i ) 

r 

i=l i<j 

Denoting overlaps by (iplaib*) = ^/yije"^'j, where yij and ipij are two real numbers, we 
have 



w 



\/\Xj\/yijyji cos [ifij - ipji 



(16) 
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For the overlap of two random states we know that the angle (p is distributed uniformly, 
while the amplitude (scaled by A^^) has an exponential distribution, p{yij) = exp {—yij). 
Therefore, for a given set of eigenvalues Aj of cta, and assuming i/ij to be independent, 
we can calculate the distribution of w using the above formula. It is a convolution of 
distributions of individual terms. In the limit r —>■ oo the central limit theorem can be 
used, resulting in a Gaussian distribution f|T3|) . 

For 10) with a small finite Schmidt rank there will be deviations from Gaussian 
because higher cumulants of p{w) ( fTTl) will be in general nonzero also in the limit 
^ oo. The above form for w (fT6!) is actually very handy for an explicit calculation 
of p{w) for 10) with a small rank r. Let us take the extreme case of rank r = 2, where 
we expect strongest deviations from a Gaussian. We therefore have only two nonzero 
eigenvalues of cta, Ai = A and A2 = 1 — A. Assuming the overlaps i/ij to be independent 
and exponentially distributed and the angles (fij to be uniform, we arrive after evaluating 
few convolutions at 



p{w) 



1 



, , Ae~~ + (1 - A)e ^ ^ H , e v'mT^ : w > 0, 



^ 4VA(l-A)+2 

The distribution is a sum of exponentials. Probability to detect entanglement, i.e., to 
measure negative value of w is V{w < 0) = 1/(4 + 2/ a/ A(1 — A)). As a function of A 
it reaches a maximum for A = 1/2, i.e., both eigenvalues of a a are equal, where it is 
equal to 1/8 (which is less than for Gaussian distribution (fT3|) ). In figure [2] we compare 
results of numerical simulation for p{w) for two cases: one with A = 1/2, for which the 
appropriate limit of eq. f|T71) gives 



p{w)x=l/2 




w > 0, 



w <0. 



and the case with an almost pure cta, A = 1/26, for which the probability of detecting 
entanglement is 5/72 ^ 0.07. Note that in the limit A — > 0, i.e., of a pure product state 
for \(f>), w will always be positive with an exponential distribution. 



3.2. Q of higher rank 

So far we have discussed only the case when Q is a one- dimensional projector, Q = 
What happens if the rank of Q is larger? If Q is of rank 2, W can be written 
as W = diWi + d2W2 with positive di^2 and Wi^2 = (|0i,2)(0i,2|)'^^- If we assume Wi 
and W2 are statistically independent, so that Wi = N"^ (iplWilip) and W2 = N^{4'\W2\4') 
are uncorrelated, then the distribution of w = diWi + ^2^2 is given by a convolution 
of distributions for Wi^2- Let us calculate the second moment of w given by the above 
sum. Using 

[ dV{ip\Wi\tlj){ilj\W2\^) = tr{WiW2) + tiWitTW2 + 0{1/N), (19) 
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Figure 2. Density of probability distribution oi w = N'^ {il;\W\ip) for random vectors 
lip) and a single W = \(f>) with having two nonzero Schmidt coefficients -s/A 
and Vl - A. We show two cases, A = 1/2 and A = 1/26, both for iV^ = 2^". Dotted 
curve is theoretical prediction ^T7\ for A = 1/26 while the dashed one is for A = 1/2, 
see eq. (|18p . We use an ensemble of 5 x 10^ random states. 



and 



we get 



tT{WiW2) = \{(pl\(p2)\', (20) 



{w-wy = dl + dl + 2did2\{(j)l\(j)2)\^- (21) 

We see that the width of the distribution of w is the same as in the case of convolution 
of two independent distributions, leading to a width df + d"^, provided the two vectors 
and 102) are orthogonal. Using similar considerations for higher moments we can 
conclude that in the case of Q = di\4>i){4'i\ having rank k, and because eigenvectors 
of Q are orthogonal, the distribution of w is Gaussian of width = Yli df, at least 
for sufficiently large k. Assuming for simplicity that all di are the same, di = 1/k, we 
obtain 

piw) = /Je-'^^--)^/^. (22) 
3.3. Mixtures of random states 

It is also interesting to consider the case in which Q is of rank 1 but we wish to detect 
the entanglement of mixed states, for instance of a mixture of m pure random states, 
as given in eq.(l2]). In such case the expectation value is w = (W^l^i) ('?/'i|), 

with W = (10) (01)"'"''. Repeating the same derivation as for Q of higher rank and pure 
states, just replacing k by m, we see that the distribution of w will be a convolution of 
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distributions for individual I?/',) and, provided these are statistically independent, the 
width of the resulting Gaussian (l22l) will be 1/m. Because the probability density p{w) 
becomes narrowly peaked about its mean W = 1 with increasing m, the probability to 
measure negative values decreases with m, namely for the Gaussian form of eq. (122|) we 
have 

„ / N 1 — erf ( \/m/2) 1 ^ /o , , 

V{w < 0) = ' ' X e""^/^ (23) 

2 v27rm 

This probability decays to zero exponentially with m. Therefore, the detection of 
entanglement for a mixture of random states is very hard. 

4. Known random state 

In this section we assume that the random state whose entanglement we want to 
measure is known in advance and furthermore, that we are able to prepare an arbitrary 
D-EW. In addition, we have to assume that our state is neither separable, nor bound 
entangled, which is true with probability which converges to one exponentially in A^. 
Therefore, for each \ip) we can prepare an optimal entanglement witness, such that its 
expectation value will be minimal. As far as D-EW are concerned, the optimal choice of 
W = Wopt is to take for Q a projector to the eigenspace corresponding to the minimal 
(negative) eigenvalue Amin of p^^, PVopt = (l^min) (^minl)"'"^- The maximal violation of 
positivity is therefore 

tr(iyoptp) = -|A^i,(pT«)|. (24) 

If we are able to measure entanglement witness with a given precisioi]|§| it is the size of 
Amin which determines the difficulty of detecting entanglement in Note that the 
optimal entanglement witness Wopt depends on the state {ip). For each state \ip) we 
have to pick a different Wopt- 

4-1. Distribution of eigenvalues of p^^ 

First, let us look at the distribution of eigenvalues after FT operation for a single 
random pure state. The eigenvalues can be written in terms of Schmidt coefficients /i, 
as A = ±fiifij. Formally, we can write the distribution of A, for A > 0, as 

dV f f dV dV 

- = ydp.ydp,5(A-p.p,)^^ (25) 

where distribution of fii can be for large N obtained from Mar cenko- Past ur law [26] for 
the distribution of r = Npf, namely 

dV/dr = v/r(4-r)/(27rr), (26) 

§ Different normalization of VF, e.g., tr = /(A^), would result in tlie maximal violation —/(A^)| Amin |. 
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by a simple change of variables t ^ fii. The result for the distribution of scaled 
eigenvalues y = NX reads 

^ = ^[(16 + l/')^(l-l/Vl6)-32E(l-yVl6)], i/G[-4,4], (27) 

where E{x) and K{x) are elliptic integrals. Note that this distribution takes into account 
only the N{N — 1) "off-diagonal" eigenvalues ±fj.ifj.j with i < j. "Diagonal" eigenvalues 
/if have the same distribution as the eigenvalues of pa and are only in number. 
The expectation value of the minimal eigenvalue equals Amin = —4/N. In fact the 
distribution of Amin becomes strongly peaked around —A/N with diminishing relative 
fluctuations as — oo. 




0.015 



Figure 3. Distribution of eigenvalues of p^^ , eq. ([2]), for mixing few states, to = 1, 2, 3 
(left plot), all for iV^ = 2^". Theoretical formula for m = 1 (^7)) is shown with the 
full curve. The right frame shows the distribution after mixing many states, m = 256, 
N"^ = 2®. Full curve is a semicircle with the center at 1/N^. All histograms are 
averages over an ensemble of 10 mixtures p 

When we mix several independent (in general non-orthogonal) random vectors, 
P — YIT the minimal eigenvalue Amin increases and the distribution becomes 

increasingly sharply peaked (for m oo we get p —>■ l/N"^ with all eigenvalues being 
equal to 1/N'^). We numerically verified the above theoretical prediction for dV/dy 
(!27|) in figure [31 We can also see that after mixing many random states (m ~ A^) the 
distribution becomes a semicircle, which is a numerical result for which we have yet no 
analytical explanation. 

4.2. Minimal eigenvalue of p^^ 

Because Amin determines the maximal violation for an optimal witness W ( HM we are 
going to look in more detail at the dependence of Amin on m and N. For a single pure 
state (m = 1) we know that the average minimal eigenvalue is Amin = — On the 
other hand, we also know that for large m, as we approach completely mixed state 
the minimal eigenvalue must scale as Amin ~ l/N"^. Therefore, the scaling of Amin has 
to change as we increase the number of mixed states m. To confirm this expectation. 
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we performed numerical simulation, calculating the average Amin for different m. The 
results are shown in figure |H Note that the average minimal eigenvalue Amin is positive 
for m > m*, with m* ^ AN'^. 
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Figure 4. Dependence of Amin on m. Overlapping of curves for different dimension N 
signals the scaling Amin ~ — 1/iV for small m in the left frame (a) and Amin ^ —l/N"^ 
for large m (but smaller than m* where Amin changes sign) in the right frame (b). We 
average over 1000,100,10, and 6 mixed random states ([2]), for N"^ — 2^,2^,2^'' and 
2^^, respectively. 



Although von Neumann entropy of a random state is large all eigenvalues of p^^ 
are very small and will therefore be hard to detect. If we assume that we are able to 
measure values of tr {W p) with accuracy e then we can, depending on the scaling of e 
with A^, tell for which values of m the detection of entanglement is possible. If e does 
not depend on i.e., precision does not increase with iV, then for sufficiently large 
iV, such that < e, detection of entanglement will be impossible. Already a single 
random state becomes from the viewpoint of entanglement detection "classical", since 
measuring a negative expectation value of its optimal entanglement witness is below 
the detection limit. If on the other hand we are able to measure e which decreases as 
1/A^, the critical merit, beyond which the entanglement detection is impossible, will be 
independent of A^, i.e., in the limit N oo the ratio merit/ N (see fig. UK). If 
however we are able to detect very small expectation values of order l/N^, then merit 
will be proportional to A^^ (see fig. |Dd). Furthermore, even with arbitrary accuracy, 
detection of entanglement with D-EW is impossible beyond m = m* oc N"^. 

Worth mentioning is that for m > N the mixed state p of eq.([2]) can not be 
used for dense coding. Indeed a quantum state is useful for dense coding [27] only if 
S{pa) — S{p) > (see e.g. [28j). For mixtures of m random states S{p) is roughly equal 
to ~ logm. On the other hand S{pa) is at most log A^. Therefore, S{pa) will be smaller 
than S{p) for m N . 
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5. Conclusions 

In this paper we have considered random states and have shown that, while their 
entanglement content is almost maximal, the detection of such entanglement is very 
difficult. This is a consequence of the complexity of a random state, which leads to a large 
number of small coefficients in the Schmidt decomposition of the state. Nevertheless, 
for random pure states, a finite success probability in the detection of entanglement 
exists also in the limit in which the Hilbert space dimension oo. This implies 

that quantum chaos alone is not sufficient to erase any trace of entanglement in the 
classical limit, provided that ideal measurements are possible. On the other hand such 
erasure becomes very efficient when coarse graining is taken into account, for instance 
when mixtures instead of pure states are considered. 

We note that all our results can be straightforwardly generalized to the case of 
unbalanced bipartition = dimT^A ^ N' = dimTie, for instance, distribution of 
w = NN' (iplWlip) is Gaussian ( fT3i) provided both dimensions N,N' are large. 

We would like to stress once more than the detection difficulties are a consequence 
of the complexity of random states. If instead one considers "regular states" such 
as the GHZ state |GHZ) = ^(|0...0) + then the optimal witness is Wopt = 

(|0min)(0min|)^^ with = ^ ( 1 0. . .0) 1 1 . . . 1) - 1 1 . . . 1) 1 0. . .0) ) which corrcspouds to the 

minimal eigenvalue Amin = —1/2 of (|GHZ)(GHZ|)'^b_ Since the value of Xmm is —1/2 
instead of —4/N as for a random state, it turns out that it will be much easier to detect 
entanglement in a "regular" rather than in a random state. This happens in spite of the 
fact that the entanglement content is larger in a random than in such a regular state. 
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